Pressure drag in linear and nonlinear quantum fluids 
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We study the flow of a weakly-interacting Bose- 
Einstein condensate around an obstacle by numerical 
solution of the Gross-Pitaevskii equation. We observe 
vortex emission and the formation of bow waves lead- 
ing to pressure drag. We compare the drag law with 
that of an ideal Bose gas, and show that interactions 
reduce the drag force. This reduction can be explained 
in terms of a 'collisional screening' of the obstacle. 

PACS numbers: 03.75.Fi, 67.40.Vs, 67.57.De 



A central issue in fluid flow concerns the origin of 
resistance or drag. In a viscous fluid, shear stresses 
induced by friction at a surface lead to skin drag. 
In an ideal fluid or superfluid, the effects of shear 
stress vanish, but normal stresses induced by pres- 
sure gradients across an obstacle still produce pres- 
sure drag. 

The recent experimental breakthroughs allowing 
the production dilute Bose-Einstein condensates 
jlj has stimulated new interest in the flow and 
modes of excitation of quantum fluids Di- 
lute Bose-Einstein condensates are compressible 
nonlinear quantum fluids, with the attractive fea- 
ture that the time-evolution can be accurately 
described by a nonlinear Schrodinger equation 
(NLSE), known as the Gross-Pitaevskii equation, 
allowing direct quantitative comparison between 
theory and experiment ||. One could envisage an 
experiment to measure the pressure drag in a di- 
lute Bose-Einstein condensate by studying the flow 
past an obstacle such as a far-detuned laser beam 
H or a foreign condensate species [||J|] . From nu- 
merical solution of the NLSE, it is known that in 
a uniform flow there is a critical velocity v c for the 
onset of drag §, however, the drag coefficient and 
exact dependence on the flow velocity remain un- 
known. 

In this paper, we address these issues by simulating 
the flow of a weakly-interacting Bose-Einstein con- 
densate around penetrable and impenetrable ob- 
jects. We show that for a flow velocity v larger 
than the critical velocity (v > v c ), the drag force 
varies quadratically with v, similar to an ideal Bose 



fluid. The principle effect of interactions is to re- 
duce drag, which we explain in terms of a 'colli- 
sional screening' of the object. Also, we show that 
the effective size of the obstacle is increased due to 
the healing length of the fluid, and consequently 
the drag is non-zero even for point-like objects. 

The NLSE in a uniform potential may be written 
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where the nonlinear coefficient, C, is determined 
by the strength of particle interactions, and V is 
the object potential. This equation is equivalent to 
two hydrodynamic equations corresponding to the 
conservation of mass and momentum, respectively, 

d t p + diJi = 0, and d t J l + djTij = , (2) 

where 



p = rmfj*ip, and J{ = —(ip*diip - ipditp*) , (3) 
are the mass density and momentum flux, and 
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is the stress tensor. The instantaneous drag force 
on an obstacle is given by 



F k (t) = / d Sj T jk (t) 



(5) 



where S defines the surface of the obstacle and 
d-Sj is an element of S in the direction of the out- 
ward normal. For a linear fluid (C = 0) and an 
impenetrable cylindrical obstacle of radius R, the 
drag law may be derived analytically: for high ve- 
locity or large object size, (v ^> fr/mR) the force 
approaches the classical limit, 



F = - Po Rv z , 



(G) 



where po is the background mass density; whereas, 
for low velocity or small objects (v <§; fi/mR), 
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F = a(mvR/h)pohv 



(7) 



where a(C) = 16/ir 2 ( 2 \H( 1} (()h[ 1} (()\ 2 is only 

weakly-dependent on velocity: is the Han- 

kel function. In this case, often encountered in 
acoustic scattering, the wave properties dominate 
and the effective object size is proportional to the 
wavelength. Inserting R ~ A = h/mv into Eq. (^), 
implies F ~ pohv, consistent with Eq. ([7]). 

For a weakly-interacting condensate (C > 0), the 
drag force must be evaluated numerically. For the 
numerical solution, we adopt the reduced units, 
t = t/ (%/tiqC), y = yj (H/^/mnoC), and tp = 
V , /v / ^o> where no is the in-flow number density. 
In reduced units, distance y, and velocity v are 
measured in terms of the healing length, | = 
h/ y/mnoC, and the speed of sound, c = \J noC/m, 
respectively; and Eq. ([!]) becomes 
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This equation is solved in 2D for both penetra- 
ble and impenetrable objects 0. A uniform flow 
in the — y-direction is imposed by multiplying the 
stationary solution by a phase gradient, e~ lv ^ . In 
principle, the instantaneous drag force can be de- 
termined by numerical integration of the stress ten- 
sor at the surface of the obstacle, Eq. However, 
for impenetrable objects, the finite grid size intro- 
duces errors in the differencing approximations to 
the surface derivatives of ip. And for penetrable 
objects, this procedure is complicated because the 
fluid-object boundary is ill-defined. The numerical 
integration can be greatly simplified by recognizing 
that the time-averaged drag force must be equal to 
the back-action on the fluid, i.e., the instantaneous 
force, 



F k (t) 
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where T defines the outer border of a simply- 
connected region of fluid, A, encircling the object. 
The second term, which corresponds to the rate of 
change of the fluid momentum within A, averages 
to zero, if the flow velocity remains constant. 

Eq. (||) may be used to calculate the drag for 
both penetrable and impenetrable objects. Fig. ([!]) 
shows a plot of the instantaneous drag, F y (t), on 
an impenetrable cylinder with radius R = 5, in a 
flow with velocity v = 1.5. Initially, the force is 
dominated by transients which depend on how the 
flow is turned on. For an instantaneous turn-on, 



reflections from the obstacle produce sound waves, 
which are subsequently absorbed at the edges of 
the box |Q. However, for longer times the time- 
averaged drag is independent of the initial condi- 
tions. Also, the time-averaged force is independent 
of the integration path T, and for barrier height, 
V > 1, only weakly dependent on the penetrability 
of the obstacle. 
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1. The instantaneous force, F y (in units of 
%*J UqC /m), on an obstacle placed in a nonlinear quan- 
tum flow with velocity v — 1.5, as a function of the 
reduced time, t. The two curves corresponds to switch- 
ing on the flow gradually (solid) or instantaneously 
(dashed). The oscillations are produced by the peri- 
odic emission of vortex pairs. Although the instan- 
taneous drag depends on the initial conditions, the 
time-averaged drag does not. 

The oscillatory behaviour of the instantaneous 
drag is produced by the periodic emission of vor- 
tex pairs. The vortex shedding frequency fol- 
lows from the phase-slip between the main flow 
and the almost stationary wake behind the ob- 
stacle: the wavefunctions for the flow and the 
wake may be written as -0 = e ~ 1 ( 1 + v / 2 )t e - lv y anc i 



respectively, where n is the mean 



density behind the obstacle (0 < h < 1, decreasing 
at higher velocity). A vortex pair is emitted each 
time the phase difference accumulates to 27r, giving 
a shedding frequency, / = (l + u 2 /2 — h)/2-K. Fig. || 
shows a comparison between the numerical results 
and the phase-slip model. The shedding frequency 
lies between the upper and lower limits set by the 
density, except at low velocity, where the shedding 
frequency falls to zero. 
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FIG. 2. The vortex shedding frequency, / (in units 
of noC/h), as a function of the flow velocity, v. The 
numerical results lie between the dashed lines, v 2 /in 
and (v 2 + 2) /Ait, which correspond to the upper and 
lower bounds predicted by considering the phase-slip 
between the main flow and an almost stationary wake. 
The error bars reflect the fluctuations in the vortex 
shedding frequency. 

The subsequent motion of the vortices is complex: 
the first pair is overtaken by subsequent pairs and 
becomes 'trapped' behind the obstacle. However, 
as apparent in Fig. 1, the vortices do not contribute 
directly to the time-averaged force, their primary 
role is to allow the separation of the wake, which 
then results in a pressure gradient across the obsta- 
cle. Fig. H shows the time-averaged drag as a func- 
tion of flow velocity. The error bars correspond to 
the residual fluctuations after averaging. The drag 
curve for the linear fluid, Eqs. (|^) and (]?]), is shown 
as a dashed line. One sees that the main effect of 
particle interactions is to reduce the drag. Below a 
critical velocity, v c , the drag in the nonlinear fluid 
falls to zero. In this region, the solutions of the 
NLSE are time independent and symmetric. The 
numeric value of v c depends on the object shape 
and penetrability. For an impenetrable cylinder, 
we find v c = (0.45 ± 0.01)c, consistent with previ- 
ous work ■ The critical velocity is larger for an 
impenetrable square barrier, and lower for a pene- 
trable Gaussian object. 
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FIG. 3. The time-averaged drag force, F v (in units 
of Ti^/rigC /rn) as a function of flow velocity, v, for a 
impenetrable cylinder with radius 7? = 3. The error 
bars indicate the magnitude of residual fluctuations 
in the time-averaged drag. The drag law for a linear 
fluid is indicated by the dashed line. The effect of the 
non-linearity is to reduce the drag which can be under- 
stood in terms of a collisional screening of the object. 
The reduced drag predicted for a screened object is in- 
dicated by the dotted line. The force as a function of 
the object size for v = 1.5 is shown inset. Note that, 
the drag is non-zero even for objects much smaller than 
the healing length. 

Above the critical velocity, a pressure drag appears 
with a velocity dependence similar to the linear 
fluid. In fact, the force can be predicted accurately 
by a semi-classical modification of the linear drag 
law, i.e., 
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F = ^PoR v 



(10) 



where R' and v' are an effective object radius and 
flow velocity, respectively. A plot of the drag force 
as a function of object size, Fig. ^( inset), indicates 
that the effective object size is extended by the 
healing length £ of the fluid. For v = 1.5, we find 
that R' = R + 0.6£. An important consequence 
of this result is that the force does not vanish for 
small objects, R <C £, e.g. a point defect in a su- 
perconductor. 



The effective flow velocity, v', may be estimated 
by considering how the flow is modified by the col- 
lisional mean-field. Above the critical velocity, in- 
coming waves are reflected by the obstacle produc- 
ing a standing wave or bow wave (see Fig. |), with 
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a maximum density variation (relative to the back- 
ground flow) of between —no and +3uq. In a semi- 
classical treatment, an oscillatory potential slows 
incoming particles by an amount corresponding to 
half the maximum barrier height, i.e., the effective 
flow velocity, v' 2 — v 2 — Ac 2 , where < A < 3. 



The dotted line in Fig. 3 is a plot of Eq. (10) us 



ing values of A obtained from the numerical so- 
lution. This 'collisional screening' model is accu- 
rate at high velocity, but less so at lower velocity, 
where the semi-classical particle treatment begins 
to break down. 

Fig. 4 show a comparison of the time-averaged 
density distribution for nonlinear (left) and lin- 
ear (right) quantum fluids. One sees that the fi- 
nite compressibility of the nonlinear fluid tends 
to suppress large density fluctuations leading to 
a smoothing of the standing wave in front of the 
obstacle. Also far downstream, the direction of 
the bow waves approaches the Mach angle, a = 
sin (c/v). For the linear fluid, c = and a = 0, 
i.e., the bow waves runs adjacent to the geomet- 
ric shadow behind the obstacle. In the nonlinear 
fluid, the 'shadow' is far less pronounced: the dark 
streak in the wake, close to the axis of symmetry, 
corresponds to the vortex street. 




FIG. 4. The time-averaged density distribution for 
v = 2.0 with (left) and without (right) interactions. 
The obstacle is an impenetrable cylinder with radius 
R — 3 centred at the origin, and the flow is from top 
to bottom. The repulsive particle interactions tend to 
smooth density variations, thereby reducing the quan- 
tum pressure experienced by the obstacle. For the non- 
linear fluid (left), the dark line in the wake, close to the 
axis of symmetry, corresponds to the vortex street. 



To summarize, we have solved the NLSE equation 
in 2D to simulate the flow of a weakly-interacting 
Bose-Einstein condensate around an obstacle. We 
observe vortex emission and the formation of bow 
waves leading to a pressure drag. We find that the 
drag force is proportional to the screened energy of 
the flow, and to the object cross-section extended 
by the effect of fluid healing. 
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